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Q>^ \ Despite encouraging advances in the design of relay codes, several important challenges remain. 

Many of the existing LDPC relay codes are tightly optimized for fixed channel conditions and not easily 
adapted without extensive re-optimization of the code. Some have high encoding complexity and some 
need long block lengths to approach capacity. This paper presents a high-performance protograph-based 
O I LDPC coding scheme for the half-duplex relay channel that addresses simultaneously several important 

issues: structured coding that permits easy design, low encoding complexity, embedded structure for 
' convenient adaptation to various channel conditions, and performance close to capacity with a reasonable 

block length. The application of the coding structure to multi-relay networks is demonstrated. Finally, 
j- ^ a simple new methodology for evaluating the end-to-end error performance of relay coding systems is 

CSJ , developed and used to highlight the performance of the proposed codes. 

o 
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L Introduction 

To harvest the cooperative gain promised by information theory [1], [2], coding schemes are needed 
that can approach the fundamental limits of the relay channel. Among the early coding results for the 
relay channel were distributed coded diversity schemes [3], [4], [5], which proposed convolutional and 
turbo codes for the fading channel under cooperation. Later on, Duman et al. [6], [7] proposed turbo 
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codes for the decode-forward (DF) relay channel in half-duplex and full-duplex modes. There is a large 
body of work focused on designing LDPC codes for the relay channel [8], [9], [10], [11], [12], [13], 
[14]. These works mostly utilize irregular LDPC codes and use density evolution (or related) techniques 
to search for optimized irregular LDPC ensembles operating at two different rates. Razaghi and Yu [11] 
produce bilayer LDPC code structures for the DF scheme. This approach has been refined in other works 
including [8], [9]. Other examples of LDPC codes for the half-duplex relay chaimel include [10], [12], 
[13], [14]. 

It has been known that the basic problem of coding for the DF relay chaimel in the power-Umited 
regime can be reduced to the following: a code and its sub-code must be designed simultaneously that 
operate with two different rates at two different SNRs at the relay and the destination, respectively. 
Despite some progress, LDPC codes for the relay channel are often painstakingly optimized to match 
to a set of channel conditions without a structure to facihtate the optimization of design, and many of 
them do not offer easy encoding. This paper proposes a class of LDPC relay codes that address three 
important issues in an integrative manner: low encoding complexity, modular structure allowing for easy 
design, and rate-compatibility so that the code can be easily matched to a variety of chaimel conditions 
without extensive re-optimization. In addition, the proposed codes offer excellent performance. 

In this paper, we concentrate on the DF protocol. For practical purposes, we further limit ourselves 
to the half-duplex relaying where the relay cannot transmit and receive at the same time. The main 
contribution of this paper is a coding scheme for the relay channel using protographs [15] that are built 
with bilayer expurgated and bilayer-lengthened structures [11]. A protograph code [15] is an LDPC code 
that can be constructed by extension from a small bipartite graph called a protograph, whose graph 
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topology is represented by a so-called proto-matrix. Protograph codes can achieve very good thresholds 
with low encoder complexity as well as fast decoding [16], [17]. The proposed approach offers flexibility 
in designing a family of rate-compatible embedded codes for relay channels. These embedded codes 
allow a coding scheme whose rate can easily be adapted to channel conditions, and are thus suitable 
for designing multi-relay coding systems. We also introduce a methodology to evaluate the end-to-end 
error performance of relay coding schemes, and demonstrate the end-to-end performance of our proposed 
codes. 

II. Background 

A. System Model 

A half-duplex single -relay channel is shown in Figure [T] Xi and Wi denote the transmitted signals 
from the source (S) and the relay (R) while Yi and Vi denote the received signals at the destination (D) 
and relay, respectively. Subscript i = 1 denotes the broadcast mode which is active for a fraction t of the 
transmission interval, and subscript i = 2 denotes the multiple access mode which is active for a fraction 
t' = I — t of the transmission interval. The received signals are: 

Vi = hsRXi + Nr^ 
Yi = hsDXi + Nd, 

Y2 = hsDX2 + hRDW2 + Nd, , (1) 

where hso and hR£i are the S-to-D and R-to-D channel coefficients, respectively, Nr-^, Nd^ and N^^ 
are the noise samples at the relay and the destination in the first and second time slots, respectively. All 
noise samples are assumed to be Gaussian with zero mean and unit variance. Ps^ = E[Xf] represents the 
source transmission power in the first time slot (duration t). Similarly Ps^ and Pr^ represent the source 
and relay transmission powers in the relayed time slot (duration t' = 1 — t). We also define SNRsr, 
SNRsD and SNRrd as signal-to-noise ratios by the relay and destination during the first time slot, 
and by the destination during the second time slot, respectively. 

B. Coding for Half-Duplex Relay Channels 

In the first time slot, a source sends a high-rate Csr^ code to a relay and a destination, using a rate that 
is decodable at the relay. In the second time slot, the relay transmits additional mutual information to help 
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the destination decode, producing an effectively lower-rate C^D^ljThus the relay coding problem consists 
of a simultaneous design of two codes that operate at two different SNRs SNRsr > SNRsd, such that 
one is a subset of the other. This problem has been attempted a number of times by forceful optimization, 
but a more streamlined approach is now available via bilayer expurgated and bilayer-lengthened LDPC 
structures [11]. In the bilayer LDPC structure, either the bit nodes or the check nodes are divided into 
two sets (called layers) and the graph, although still a bipartite graph, is re-drawn with two rows of check 
nodes sandwiching the bit nodes, or vice versa [11, Fig. 3 and 6]. This provides a convenient way to 
build and illustrate subcodes, which we use for our purposes. For a more comprehensive treatment and 
justification of the bilayer structure please see [11], [13]. 



C. Protograph Codes 

A protograph [15] is a Tanner graph with a relatively small number of nodes. A protograph code is 
an LDPC code constructed from a protograph by a copy-and-permute operation, where the protograph 
is copied N times and the N variable-to-check pairs (edges) corresponding to the same edge type of 
the original protograph are permuted. Thus the protograph code has the same rate and the same node 
distribution as the original protograph. Several capacity -approaching protograph code designs have been 
proposed in [17], [18]. Figure |2] shows a rate- 1/2 protograph from [19] consisting of 7 variable nodes and 
4 check nodes that are interconnected by 24 edge types. The dark circles represent transmitted variable 
nodes, the white circle is a punctured node and the circles with a plus sign are parity check nodes. This 
protograph can be represented by the proto-matrix 



H 



1/2 



(2) 



'^l 2 1 0^ 
3 11110 
1 2 2 2 1 1 
y0 200002y 

where the rows and columns represent the check nodes and variable nodes respectively. The entries of the 
matrix represent the number of parallel edges that connect the variable node and the check node. Using 
the PEXIT technique [20], the threshold of this code is 0.439 dB, a gap of 0.252 dB of capacity. The 



In our formulation the source may also be active in the second time slot, but it will be constrained to synchronously transmit 
the same signal/codebook as the relay. Thus in this work the distinction between orthogonal and non-orthogonal transmission 
is only in the received SNR during relay time slots. Since we parameterize our analysis by received SNR without any direct 
reference to source and relay powers, both orthogonal and (synchronous) non-orthogonal transmissions are subsumed in the 
following developments. 
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linear minimum distance of this code grows linearly with codeleng tH3 which is necessary for avoiding 
an error floor [21]. 

III. Design of bilayer protograph-based LDPC codes for the relay channel 

A. Bilayer Lengthened Structure 

A bilayer lengthened structure for the relay channel [11] is shown in Fig. [T] where a high-rate code is 
constructed by adding variable nodes to the graph of a low-rate code. Geometrically, the overall graph 
contains two layers (sets) of variable nodes and one set of check nodes. One layer of variables connecting 
to the checks forms a graph of a lower-rate code, and the entire graph constitutes a high-rate code. The 
corresponding parity check matrix has the form: 

HsR, = [HsdM , (3) 

where HsDi is the parity matrix of the lower-rate code which is equivalent to the sub-code (Layer 1 or 
CsDi) and H^. is the extension matrix representing Layer 2. Layers 1 and 2 combined create a capacity- 
approaching code for the source-relay link {CsrJ as shown in Fig. [T] The relay decodes this codeword, 
then protects the codeword bits in Layer 2 via k2 parity bits. These parities are encoded with Crd^ and 
transmitted to the destination during the second time slot. The destination, using these k2 parity bits, can 
reliably detect the bits in Layer 2, therefore they can be eliminated from the code graph. The remaining 
graph contains just Layer 1, which constitutes a capacity approaching code at the SNR of the SD link, 
which is decodable at the destination. 

We will describe in the following an example for the design of the bilayer lengthened codes with rates 
R = n = 1, 2. . . ., as an extension of the rate-1/2 protograph of Eq. For each new code, three 

^If a random permutation per each edge of its protograph is used. 
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variable nodes (or columns) are added, as follows: 
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where Xi and variables denote the number of parallel edges in the extension graph, to be determined. 
Variables yi in the first row correspond to the check node that connects to the degree- 1 variable node. 
Variables xi in rows 2-4 correspond to the remaining check nodes. In order to preserve the linear growth 
of minimum distance for the bilayer codes, the column sums in rows 2-4, namely sums of columns 
designated with variables Xi, should be 3 or higher [21]. 

We further simplify the problem by setting the maximum number of parallel edges to 2 (i.e. Xj, y^ = 
0, 1, 2). These constraints limit the range of parameters thus simphfying the optimization. Our cost 
function is the threshold which is calculated via the PEXIT method. For this example, the best threshold 
is given by: 
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TABLE I 

Thresholds of proposed protograph codes 



Code 


Protograph 


Capacity 


Gap to 


Rate 


threshold (dB) 


threshold (dB) 


capacity 


Bilayer lengthened design 


1/2 


0.439 


0.187 


0.252 


2/3 


1.223 


1.059 


0.164 


3/4 


1.720 


1.626 


0.094 


4/5 


2.136 


2.040 


0.096 


5/6 


2.455 


2.362 


0.093 


6/7 


2.718 


2.625 


0.093 


7/8 


2.941 


2.845 


0.099 


8/9 


3.125 


3.042 


0.083 


9/10 


3.295 


3.199 


0.096 


Bilayer expurgated design 


3/4 


1.720 


1.626 


0.094 


2/3 


1.182 


1.059 


0.123 


7/12 


0.809 


0.590 


0.219 


1/2 


0.420 


0.187 


0.233 


5/12 


0.144 


-0.185 


0.329 


1/3 


-0.263 


-0.497 


0.234 



The iterative decoding thresholds of these codes calculated by PEXIT technique [20] are given in 
Table in For rates > 2/3, the produced codes have iterative decoding thresholds within 0.1 dB of the 
capacity, and the rate-2/3 code has a threshold within 0.152 dB. 



B. Design of Bilayer Expurgated Protograph Codes 

A bilayer expurgated structure for the relay channel [11] is shown in Fig. [3] where a low -rate code 
is constructed by adding check nodes to the graph of a high-rate code. The bilayer graph in this case 
contains two layers (sets) of check nodes and one set of variable nodes (unlike the bilayer lengthened 
structure). The variable nodes together with one layer of check nodes form the graph of a high-rate code, 
and the entire graph constitutes a low-rate code. The corresponding parity check matrix has the form 

HsRi 



HsDi 



(5) 
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Layer 2 




Fig. 3. Bilayer expurgated protograph design for half-duplex relay with Rsdi ~ and Rsr^ ~ , N is the 

number of variable nodes 

where Hsr^ is the parity check matrix of a high-rate LDPC capacity-approaching code for the source- 
relay Unk (representing the sub-code CsRi of Layer 1), and He is the extension matrix representing 
Layer 2. Layers 1 and 2 together create a capacity-approaching code CsDi for the source-destination 
Unk as shown in Fig. |3] The source transmits a CsBi codeword. The relay, after decoding the source 
codeword, produces k2 additional "parity" bits using the extension matrix H^, re-encodes these k2 bits 
with a codebook Crd2 ^^d transmits to the destination. At the destination, the ^2 parity bits that are 
reliably detected essentially provide k2 additional check values on the source codeword from Csr^ ■ This 
is equivalent to decoding a Csdi codeword (with a lower rate) at the SNR of the source-destination link. 

We now describe an example for the design of bilayer expurgated codes. In the interest of brevity, 
we only present one construction of bilayer codes. However, the proposed method is completely general; 
we can start from any high-rate code to build a family of expurgated codes. The starting point for 
the following construction is the rate-3/4 protograph code designed in the last subsection. The rate-3/4 
protograph contains 13 variable nodes (one of them punctured) and 4 check nodes. Low-rate bilayer 
codes expurgated from this code have rates R = ^^^rT^> where n is the number of checks added. For 
each value of n, the new code will require a search. As a representative sample, we concentrate on the 
search for n = 1, which yields a rate-2/3 code. The search space for this new code is in the form 
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H. 



search 
2/3 



3/4 



(6) 



\ y xi X2 . . . xu J 
where y variable denotes the number of parallel edges connecting to the punctured variable node and 
Xi variables denote the number of parallel edges connecting to remaining variable nodes except the 
degree- 1 variable node. Having a degree- 1 variable node improves the iterative decoding threshold of the 
protograph code [17]|^ 

The optimization process is simplified by further limiting the maximum number of parallel edges. For 
example, we can limit y G {1,2} and Xj G {0, 1,2}. Proto-matrices of new expurgated protographs are 
given by 



H. 



2/3 



H. 



3/4 



Y0100000110012y 



-^7/12 



H: 



2/3 



y0100100100002y 



^1/2 



H. 



7/12 



\o2ooooooioiooy 



It is well known that degree- 1 nodes are not to be used in irregular LDPC codes, but their effectiveness in structured LDPC 



codes is well documented [22]. 
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h: 



5/12 



^1/3 



H 



1/2 



^0200000011000/ 
/ \ 



■"7/12 



\02ooiooooooooy 

The iterative decoding thresholds of these codes are given in Table IH The produced codes all have 
decoding thresholds within 0.23 dB of the capacity except the rate-7/12 one which has a decoding 
threshold within 0.33 dB of the capacity. 



IV. Numerical results 

A. Coding thresholds 

We design nested bilayer lengthened and expurgated protograph codes for the half-duplex relay channel, 
which can operate at any two rates with iterative decoding thresholds close to capacity, as shown in 
Table H For comparison, the thresholds of bilayer lengthened and bilayer expurgated irregular LDPC 
codes proposed in the literature are shown in Tabled All codes use a rate pair {RsDi — 0.5, RsBi — 0.7). 
Our proposed bilayer lengthened and bilayer expurgated codes are better than the codes proposed in [14] 
and within 0.1 dB of codes reported in other worksQ Thus, a simple design approach has yielded a 
family of codes that offer a performance comparable with other highly-optimized bilayer LDPC codes, 
while offering important advantages in terms of rate-compatibility, low encoding complexity, and the fast 
decoding properties of protograph codes. 



B. Simulation Results 

So far we have represented codes only in the form of proto-matrices (protographs). As mentioned 
in Section III-C[ a protograph code (an equivalent LDPC code) is a large derived graph constructed by 
copy-and-permutation operation on protograph, a process known as lifting. 



''The rate of [11], [14] is slightly less than 1/2, thus slightly skewing the comparison in their favor. 



February 20, 2013 



DRAFT 



11 



TABLE II 

Comparison of thresholds (Ri ~ 0.5, R2 ^ 0.7) 





Lengthened 


Expurgated 


[11] 


[14] 


[8] 


Our BL code 


[11] 


[9] 


[8] 


[14] 


Our BE code 


Gapi 


0.164 


0.3854 


0.1039 


0.239 


0.514 


0.258 


0.284 


0.6323 


0.233 


Gap2 


0.120 


0.1758 


0.0945 


0.152 


0.084 


0.084 


0.084 


0.215 


0.123 



Our protograph codes are derived from protographs in two lifting steps. First, the protograph is lifted 
by a factor of 4 using the progressive edge growth (PEG) algorithm [23] in order to remove all parallel 
edges. Then, a second lifting using the PEG algorithm was performed to determine a circulant permutation 
of each edge class that would yield the desired code block length. 

In our nested bilayer lengthened codes of Section ITlI- A [ the parity check matrix for the lower-rate code 
can be obtained by removing certain columns from the parity check matrix of the higher rate code, and 
this produces economies in the design of the decoders. In fact, it is enough to design a decoder for the 
largest rate code (9/10). To decode the lower-rate codes, the missing coded bits are replaced by erasures 
at the decoder. In the same manner, the bilayer expurgated codes generated in Section IIII-BI only need 
the decoder of the lowest rates (1/3). Other higher-rate codes are decoded by ignoring redundant checks 
at the common decoder. 

The performances of our bilayer lengthened protograph codes with rates 1/2, 2/3 and 3/4 over a 
binary-input additive white Gaussian noise (BI-AWGN) are shown in Fig. |4l In addition, performances 
of bilayer expurgated protographs with rates 1/3, 5/12, 1/2, 2/3 and 3/4 are shown in Fig. |5] All codes 
are simulated with the information block-length of 16k. To construct codes with this 16k information 
block-length, the protographs for the lengthened codes are first lifted by a factor of 4, and subsequently 
again by factors 1365, 683, and 455, for codes of rate 1/2, 2/3, and 3/4, respectively. A similar two-step 
lifting is used for the expurgated codes. 

The decoder is a standard message passing decoder where the maximum number of iterations is set 
to 200. Log-likelihood ratio (LLR) clipping and other decoding parameters are according to [24]. No 
error floors were observed down to the word error rate (WER) 10~^. The gap to capacity of these 
codes is within 0.6 dB of their iterative decoding threshold and within 0.8 dB of their capacity limits at 
WER « 3 X 10-*^. 

We now consider a half-duplex relay channel with information block-length of 16380 bits. We assume 
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0.5 1 1.5 2 2.5 3 

SNR (dB) 



Fig. 4. Performance of rate 1/2, 2/3, 3/4 bilayer lengthened protograph codes designed in Section UlI-AI with information 
blocklength of 16k 

the SR link supports rate Rsr^ =3/4 and the SD link supports rate RsDi = 1/2- For this example, 
the time division t = 0.75 is chosen for the source and relay. The rate needed in the RD link is 
Rb.D2 =3/4 [13, Eq. 5]. This is the same as the rate of SR code but with a different codeword-length. 
Under these conditions, the achievable rate of the relay channel is 0.5625 [13, Eq. 3]. 

We now consider the bilayer expurgated coding scheme (Fig. |3]l with information blocklength 16380. 
CsRi is implemented by the rate-3/4 protograph given in Section ITlI-AI The code Crd^^ which protects 
additional parities generated at the relay, has information blocklength 5460. Crd2 is constructed from 
the same protograph as Csr^ , but with a shorter length. The destination decodes a Csr^ codeword plus 
these additional parity values, which constitutes a codeword of CsDx in the same manner as [11]. The 
word error rate (WER) and bit error rate (BER) performances of these codes are shown in Fig. |6] The 
shorter block length of the rate-3/4 RD code has resulted in a wider waterfall region, highlighting a 
weakness of the successive decoding approach at smaller block lengths. 

The bilayer lengthened structure involves four codes (Fig. [T]). Using the above relay channel parameters, 
CsRi has the information block-length 16380 and codeword block length 21840. Each codeword of Csr^ 
is composed of a CsDx codeword with blocklength 10920 and an extension of 10920 bits derived via 
the extension matrix from Eq. Q. The relay, after decoding and recovering the codeword, multiplies 
the extension bits by the parity check matrix of a rate- 1/2 code Ci to calculate a syndrome of length 
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Fig. 5. Performance of rate 1/3, 5/12, 1/2, 7/12, 2/3, 3/4 bilayer expurgated protograph codes (Section IlII-BK information 
block-length 16k 




0.5 1 1.5 2 2.5 

SNR (dB) 



Fig. 6. Performance of component codes for a relay channel: i?sDi=l/2, RsR^^il^ and Rrd^ ~ 3/4 using bilayer expurgated 
structure (Section HIFbI i with a = 1.4 dB and /3 = 1.6 dB 
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0.5 1 1.5 2 2.5 3 

SNR (dB) 



Fig. 7. Performance of component codes for a relay channel: i?sDi = l/2, Rsri=3/4, Rc^ = 1/2 and Rrd2 = 3/4 using 
bilayer lengthened structure in Section UlI- Al with a = 1.4 dB and /? = 1.6 dB 

5460. Ci is lifted from the rate- 1/2 protograph given in Eq. The syndrome is then encoded by Crd2' 
transmitted and decoded at the destination to recover the syndrome. Csr^ and Crd^ are similar to their 
counterparts in the bilayer expurgated scheme. Also, CsDi and Ci have the same code structure and 
blocklength. The WER and BER performances are shown in Fig. |7] Performance of CsDx in the bilayer- 
lengthening structure is worse than in the bilayer expurgated structure due to its shorter blocklength. 

Although the example above was coached in terms of orthogonal transmissions, it also applies to a 
non-orthogonal system with correspondingly lower source and relay powers (see footnote [T}. 

C. End-to-End Performance of Relay Coding Systems 

Many works in contemporary literature represent the performance of relay channel codes by illustrating 
two component codes: the source-relay code, and the source-destination code. In this approach it is 
implicitly assumed the relay-destination code is an ideal code. However, the relay-destination code often 
operates at a smaller block length, thus it has a wider waterfall and may very well be the bottle neck 
for the entire system, therefore providing the thresholds or even simulations for only the two other sub- 
codes may not always be fully illuminating of the overall performance. We believe there is need for a 
comprehensive end-to-end performance metric. 

In general the end-to-end relay channel error is a function of the three SNRs of its constituent channels. 
But showing this dependence requires a four-dimensional plot, which is not practical. We construct a 
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simple but useful error plot by assuming 

SNRsR = SNRsD + « SNRuD = SNRsd + /3 , (7) 

where a and /3 are fixed constants. The affine relationship of the SNRs allows them to be displayed on the 
same axis, producing a simple end-to-end error plot which can be thought of as a 2-dimensional slice of 
the general four-dimensional error characterization of the relay channel, a and (3 can be carefully chosen 
(as shown below) to generate a reasonable and insightful overall representation of the performance of 
the relay channel that does not ignore or trivialize any of the parameters of the system. 

To simplify the simulations, we develop a tight bound on the end-to-end error based on component 
errors in the system. For the purposes of exposition we concentrate on the expurgated-bilayer codes 
(Figure |3]l, where three types of errors may happen in the system: Er is the error event at the relay, Ejn) 
is the error in the decoding of the "extra parity" arriving from the relay to the destination, and finally, 
Ed is the error event at the destination in the final decoding of the source message, and the complement 
of event E is shown with E. The bound is shown in Eq. (fTOl ). 

Pe = P{Ed\E^)P{Er) + P{Ed\Er)P{Er) (8) 

= [P{Ed\E^,E^)P(E^) + P{ED\EnD,E^)P{ERD)]P(E^) + P{Ed\Er)P{Er) (9) 

< P{Ed\E^,'Er) + P{Erd) + P{Er) (10) 

It is easy to see the bound is tight because component codes are used in a regime where their errors 
are, conservatively, no more than 10^^ ~ 10^'^, therefore with a very good approximation P{Er) f» 1, 
P{E^) « 1, P{Ed\Er) ^ 1 and P{Ed\Erd,Er) ^ 1. 

We are now ready to calculate the values of a and /3, which will determine how the individual code error 
characteristics are combined to produce the end-to-end performance. We are interested in cases where 
none of the component errors dominate the others, because when one of the link errors dominates, the 
performance will be essentially determined by the code on the dominating link, which is already known 
from the literature on point-to-point channels. So the interesting case is when no link errors dominate, 
i.e. the waterfall regions of the component curves coincide at their starting points (approximately the 
threshold). 

The performance of two bilayer codes with a rate pair (Rsd = 1/2, RsR. = 3/4), a = 1.4 dB, /3 = 1.6 
dB and time sharing t = 0.75 are shown in Fig. [8] To understand how far is the end-to-end performance 
of our codes from theoretical limits, we follow the convention of point-to-point channels and produce 
a value of SNR at which a random code of infinite length and rate similar to the code under study is 
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Fig. 8. End-to-end performance bounds for relay coding schemes in Fig.[7]and Fig.|6]with a = 1.4 dB and P = 1.6 dB 

capable of theoretically supporting error-free communication. Similarly to the point-to-point case, this 
involves inverting the capacity formula by inserting the rate and extracting the corresponding SNR. To 
produce a single SNR value, we use Eq. d?]), assume that each of the component codes is at the rate that 
the mutual information of the link supports, and use the relay mutual information formulas in [13]. Using 
these parameters, the relay channel of the above example has a throughput of 0.5625 and the code rates 
correspond to the limiting SNR of SNRsd = 0.225 dB. As seen in this figure, the gap-to-capacity of 
bilayer expurgated and bilayer lengthened relay coding schemes are about 0.7 dB and 1.2 dB respectively 
at the frame error rate (PER) of 2 x 10~^. The end-to-end error bound of the bilayer lengthened code 
is worse than that of the bilayer expurgated code because of the block-length issues. This phenomenon 
has also been reported by several previous works including [11], [14]. 



In a practical two-relay scenario (Fig. |9]l it is likely that the relays do not have precisely and deter- 
ministically identical channels to the source, destination, and each other. Therefore one relay is likely to 
be "stronger" and decode first, then this relay will be able to help another relay to decode, and then the 
two relays together will assist the destination. 

We consider transmissions in three time slots for the source and two relays. The transmitted signals 
from the source (S) and two relays (Ri and R2) are denoted with X, W and Q, respectively, and the 
received signals at the destination and two relays are denoted with Yi, Vi and Pi, respectively, where i 
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Fig. 9. The expurgated coding structure for tlie half-duplex two-relay channel 

indicates the time slot index. In the first time slot, the source encodes its message using the code Csr^ ■ 
The first relay decodes the source message, but the second relay and the destination cannot (yet). During 
the second time slot, the first relay generates k2 parity bits using the sub-graph denoted Layer l]^ then 
encodes these parity bits using another LDPC code Ci. The second relay, after decoding the /c2 parity 
bits, decodes the source message. Then, the second relay computes parity bits using the subgraph 
denoted Layer 3, encodes the bits with another LDPC code C2, and transmits to destination. The 
destination decodes Ci, C2, and finally the source message with the help of k2 + k^ additional parity 
bits from two relays. The achievable rate using the above strategy is a special case of the achievable rate 
in [25]. 

For demonstration, consider an example where Ri, R2 and D can reliably receive source signals 
with a rate Rsr^ = 3/4, Rsr2 = 7/12 and Rsd = 1/3, respectively. The Ri-to-D channel supports 
RriD = 1/2. Assuming ti = 0.6, ^2 = 0.2 and ts = 0.2, the achievable rate of the two-relay channel is 
0.45 [13]. The rate of Ci and C2 are 1/2 and 3/4, respectively. Information block length is 16380. Csr^, 
CsR2 and Csd are constructed from protographs given in Section HII-BI Csr2 is rate- 7/ 12 and Csd is 
rate- 1/3, resulting in k2 = 3640 and k^ = 5460. Although the explanation of this example was in terms 
of orthogonal transmissions, the precise same example also applies to a non-orthogonal (beamforming) 
relay channel with correspondingly lower source and relay power. 

^Like the one-relay case, these bits are the syndrome of a parity check extension matrix. 
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SNR (dB) 

Fig. 10. Perfoimances of component codes used in a two-relay channel: Rsr^ = 3/4, Rsr2 = 7/12, Rsd = 1/3, Rcx = 1/2 
and i?c2 = 3/4. Csrx, Csr-^ and Csd all have codeword blocklengths of 23660, while Ci and C2 have codeword blocklengths 
of 7280. 

C\ arises from the protograph in Eq. (|2) and C2 is the same code as Csr^, but with a shorter 
blocklength. Fig. [TOl shows the performance of CsBi^ Csr-i^ Csd, Ci and C2, which operate within 0.6 
dB, 0.7 dB, 0.9 dB, 1.5 dB and 1 dB, respectively of their capacity Umits at WER = 10^^. As expected, 
Ci is the bottleneck in this two-relay channel because of its short blocklength. We note that this gap is 
an outcome of a relatively short overall block length and a successive decoding approach, both of which 
are practical considerations and can be relaxed. 

Remark 1: In the multi-relay scenario, there are multiple options for the decoding and transmission at 
relays. For example, two relays can cooperatively beamform towards a destination. The required codes in 
this case are fundamentally similar to the one-relay channel, therefore we do not consider them separately 
in this paper 

Remark 2: In the multi-relay scenario the transmission time is divided multiple times which shortens 
the component block-lengths. Therefore we did not pursue a multi-relay generalization of the bilayer 
lengthened structure due to its susceptibility to short block length effects. 

Remark 3: Optimization of the time division to maximize the achievable rates has been pursued in 
papers on the capacity of relay channel [26] and can be extended to the multi-relay scenario. Depending 
on the channel gains, the optimal length for one of the time slots may be zero, in which case one of the 
relays must be shut off. 
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Remark 4: The developments throughout this paper apply to either orthogonal relaying or to non- 
orthogonal relaying where concurrent transmissions use the same codebook. 

VI. Discussion and Conclusion 

This paper presents a simple approach for constructing relay coding schemes based on bilayer length- 
ened and bilayer expurgated protograph codes which perform within a fraction of dB of the capacity. 
The proposed codes allow easy design, flexibiUty in matching to various relay chaimel conditions and 
low encoding complexity and can be extended to multi-relay networks. A framework for end-to-end 
performance evaluation of the relay codes is also provided. 

Nested protograph codes have also been used for the point-to-point channel [27] where, like the 
present work, multi-component protograph codes have been considered. However, in relay channels the 
component codes are emitted by different transmitters, unlike the point-to-point case, and see different 
channels with different SNRs. Thus, relay networks represent a different scenario whose challenges are 
addressed by a bilayer structure [11] for the overall code and a coset code at the relay, techniques that 
have no counterpart in the point-to-point case. 
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